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ABSTRACT 

Light-front analysis of vr” mesons in Mg-Mg collisions is carried out. The phase space 
of secondary pions is naturally divided into two parts in one of which the thermal equilibra¬ 
tion assumption seems to be in a good agreement with data. Gorresponding temperatures 
are extracted and compared to the results of other experiments. The experimental results 
have been compared with the predictions of the Quark Gluon String Model (QGSM) and 
satisfactory agreement between the experimental data and the model has been found. 
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Keywords: NUCLEAR REACTION Mg(Mg, tt" X), at P=4.3 A GeV/c; 

measured pion distributions; the analysis in light front variables 
is carried out; deduced thermal equilibration, which is characterized 
by the temperature T; Gomparison with Quark Gluon String Model. 
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1. INTRODUCTION 


In the experiments with beams of relativistic heavy ions one hopes to observe the 
extreme conditions when the phase transitions in nuclear matter are expected (see, e.g. 
[1-3]). 

For the experimental study of such transitions it is necessary to understand the mecha¬ 
nism of collisions and investigate the characteristics of multiparticle production in nucleus- 
nucleus interactions. The study of single particle inclusive processes is one of the simplest 
and effective tools for the understanding of dynamics of multiple production of secon¬ 
daries. 

In this respect it is important to investigate the properties of 7r“ mesons, which are 
predominantly produced particles carrying the information about the dynamics of collision 
and which are reliably identihed. Besides, the pion production is the predominant process 
at Dubna energies. 

Our previous results on pion production experiment (cross-sections, multiplicities, ra¬ 
pidities, transverse momenta, intercorrelations between various characteristics, etc) using 
the streamer chamber spectrometer SKM-200 and its modified version GIBS in inelastic 
and central nucleus-nucleus interactions are presented in [4-6]. 

In this paper we present the light front analysis of tt” mesons produced in Mg-Mg col¬ 
lisions. It is assumed that interactions of identical nuclei ( Mg-Mg) give the possibility of 
better manifestation of nuclear effects than the interactions of asymmetric pairs. In some 
cases light front analysis [7] seems to be more sensitive to the details of the interaction 
mechanism as compared to the presentation of data in terms of the well known Feyman 
variables xp, rapidity Y etc. 


2. EXPERIMENT 


The data were obtained using the SKM-200 facility and its modified version GIBS of 
the Dubna Joint Institute for Nuclear Research. SKM-200-GIBS consists of a 2m streamer 
chamber, placed in a magnetic held of ~ 0.8 T and a triggering system. The streamer 
chamber was exposed by beam of Mg nuclei accelerated in the synchrophasotron up to a 
momentum of 4.3 GeV/c per incident nucleon . The solid target (Mg) in the form of thin 
disc with thickness 1.56 g/cm^ was mounted within the hducial volume of the chamber. 
The triggering system allowed the selection of ’’inelastic” and ’’central” collisions. 
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The inelastic trigger selected all inelastic interactions of the incident nuclei on the 
target. The central trigger selected events dehned as those without charged and neutral 
projectile spectator fragments ( P/Z > 3 GeV/c) within a cone of half angle 0c/i=0n=2.4° 
(the trigger efficiency was 99% for events with a single charged particle and 80% for events 
with a single neutron in the cone). The trigger mode for each exposure is dehned as 
T (0c/i ,0n ) (0c/i and 0^ expressed in degrees and rounded to the closest integer value). 
Thus Mg-Mg interactions obtained on the set-up correspond to the trigger T(2,2). The 
fraction of such events was ~ 4 ■ 10“^ among all inelastic interactions. The experimental 
setup and the logic of the triggering systems are presented in Fig.l. 

Primary results of scanning and measurements were biased due to several experimental 
effects and appropriate corrections were introduced. The biases and correction procedures 
were discussed in detail in [4,6]. Average measurement errors of the momentum and 
production angle determination for tt” mesons are < AP/P >= 1.5%, A0 =0.1°. 


3. LIGHT FRONT PRESENTATION OF INCLUSIVE DISTRIBUTIONS 


An important role in establishing of many properties of multiple production is played 
by the choice of kinematical variables in terms of which observable quantities are presented 
(see in this connection, e.g. [8]). 

Here we propose unihed scale invariant variables for the presentation of single particle 
inclusive distributions, the properties of which are described below. 

Consider an arbitrary 4-momentum and introduce the light front combina¬ 

tions [9]: 


P± = Po ± Ps 


( 1 ) 


If the 4-momentum p^ is on the mass shell {p^ = m?), the combinations pj_, (where 
Pt = (Pi) P 2 )) dehne the so called horospherical coordinate system (see, e.g. [10]) on the 
corresponding mass shell hyperboloid Po — P ^ 

Let us construct the scale invariant variables [7]: 




± 


P± 

p± +p± 


( 2 ) 


in terms of the 4-momenta p)) of particles a, b, c, entering the inclusive reaction 

a + b ^ c + X. The 2 ;-axis is taken to be the collision axis, i.e. Pz = Pa = Pl- Particles 
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a and b can be hadrons, heavy ions, leptons. The light front variables in the centre of 
mass frame are dehned as follows [7]: 




. E±p^ _ E + \p^ 

Vi “ Vi 


( 3 ) 


where s is the usual Mandelstam variable, E = \JpI+ Pt E and pz are the energy 
and the - component of the momentum of produced particle. The upper sign in Eq. 
(3) is used for the right hand side hemisphere and the lower sign for the left hand side 
hemisphere. It is convenient also to introduce the variables 


= Tlnle^ 


in order to enlarge the scale in the region of small 

The invariant differential cross section in terms of these variables looks as follows: 

^ ^ da ^ I da 

dp TT d^^dp\ TT d(^dp'^ 

In the limits of high pz {\pz\ 3> pr) and high pt {pt \'Pz\) the variables go over to the 
well known variables xp = ‘Ipzj^fs and xp = 2pt/\/s, respectively, which are intensively 
used in high energy physics, ,^^-variables are related to Xp, Xp and rapidity y as follows: 

\ ± \/xl + x^^ ; Xp = (5) 

y = ±7^ In- 2 — ]mT = \IPp + m^ (6) 

The principal differences of distributions as compared to the corresponding xp ~ 
distributions are the following: 1) existence of some forbidden region around the point 
= 0, 2) existence of maxima at some in the region of relatively small |.^^|, 3) 
existence of the limits for |,^^| < m/y/s. The maximum at is also observed in the 
invariant differential cross section — -r^. However, the region 17^1 > 17=*" I goes over to 

TT dC^ 

the region |(^=*=| < |(^=*=| and vice versa. 

Note that the light front variables have been introduced long time ago by Dirac [9] 
and they are widely used now in the treatment of many theoretical problems (see, e.g. 
original and review papers [11-19] and references therein). They have been used also in a 
number of phenomenological applications (see, e.g. [20]). 
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4. THE ANALYSIS OF PION DISTRIBUTIONS IN TERMS OF LIGHT 

FRONT VARIABLES. 


The analysis has been carried out for the 7r~ mesons from central (trigger T(2,2)) 
Mg-Mg collisions (~ 6200 events, ~ 50 000 7r“ mesons). In Figs.2 and 3 xp - and 

~ distributions of all 7r~ mesons are presented. The experimental data for invariant 
distribution (I/tt) ■ dN/d(^ are shown in Fig.4. The curve is the result of the polinomial 
approximation of the experimental distribution and the maximum is observed at = 
2.0 ± 0.3. The value is the boundary of the two regions with significantly different 
characteristics of secondaries. 

In Figs.5 and 6 the pp and the angular distributions of 7r~ mesons in the forward 
hemisphere in different regions of C"*" ( C"*" > C’*' ^md C"*" < C’*') ^ire presented. Similar 
results have been obtained for the backward emitting tt” mesons. 

One can see from Figs.5 and 6, that the p\ and angular distributions of tt” mesons 
differ signihcantly in < ,^+) and < C’*' (■C"'" > ■C"'") regions. The angular 

distribution of pions in the region (C’*' < (C+ is sharply anisotropic in contrast to the almost 
flat distribution in the region > (C+. The slopes of pp ~ distributions differ greatly 
in different regions of C'*'. The average values < p\ > in these two regions also differ: 
< p^ >= (0.027 ± 0.002) (GeV/c)^ in the region > (+; < p^ >= (0.103 ± 0.009) 
(GeV/c)2 in C+ < C+- 

The flat behaviour of the angular distribution allows one to think that one observes a 
partial thermal equilibration in the region |,^^| < > |(C=*=|) of phase space. 

Note, that the paraboloids of constant 


separates two groups of particles with signihcantly different characteristics. 

To describe the spectra in the region > (C+ the Boltzmann f{E) ~ and 

the Bose-Einstein (B-E) /(F) ~ (e“^/^ — 1)“^ distributions have been used. 

The distributions (I/tt) ■ dNjdCp^^ dN/dp^, dN/dcosQ look in this region as follows : 


da 

dpp 

da 

dcosO 
1 da 



fi.E)dp, 


rPmax , ^ n 

/ f{E)p^dp 
Jo 

[Pt. max „ „ / , o 

/ Ef{E)dpl 
Jo 


( 8 ) 

( 9 ) 

( 10 ) 
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where: 


PT,^ax = - ml 

Pz,max = [Pt + m^ - (^+v^)2]/(-2^+V^) 

Pmax = {-^^^/scosQ + \J (^+v^)2 - mlsiv?Q)/siu^Q 

The experimental distributions in the region > (■+ have been htted by the expres¬ 
sions (8), (9), (10), respectively. The results of the £t are given in Table 1 and Figs. 5 -i- 
7 (by the Boltzmann distribution). They show a rather good agreement with experiment. 
One can see from the Table 1, that the values of the parameter T extracted by htting the 
data with Boltzmann and Bose-Einstein distributions coincide within the errors. 

In the region C,'^ < the ~ distributions has been fitted by the formula 

~ a ■ -h (1 - a) • ( 11 ) 

and the ^ distributions by the formula 

r ^ ~ 

The results of the fit are given in Table 2 and Figs.5 and 7. 

Thus the spectra of 7r“ mesons in the region > (C+ are satisfactorily described 
by the formulae which follow from the thermal equilibration. The same formulae when 
extrapolated to the region C,'^ < (^+ (Fig.7) deviate significantly from the data. On the 
other hand the dependence (1 — .^^)"' is in a good agreement with experiment in the region 
and deviates from it in the region > (C+ (Fig.7). 

Several theoretical models of nucleus-nucleus collisions at high energy have been pro¬ 
posed [21]. In this paper Quark Gluon String Model (QGSM) [22] is used for a comparison 
with experimental data. The QGSM is based on the Regge and string phenomenology of 
particle production in inelastic binary hadron collisions [23]. The QGSM simplifies the 
nuclear effects (neglects the potential interactions between hadrons, coalescence of nucle¬ 
ons and etc.). A detailed description and comparison of the QGSM with experimental 
data in a wide energy range can be found in paper [24]. To describe the evolution of 
the hadron and quark-gluon phases, this model uses a coupled system of Boltzmann-like 
kinetic equations. Nuclear collisions are treated as a mixture of independent interactions 
of the projectile and target nucleons, stable hadrons and short lived resonances. QGSM 
includes low mass vector mesons and baryons with spin 3/2, mostly A(3/2,3/2) reso¬ 
nances. The procedure of generation consists of 3 steps: the definition of configuration of 
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colliding nuclei, production of quark-gluon strings and fragmentation of strings (breakup) 
into observed hadrons. The formation time of hadrons is also included in the model. The 
QGSM has been extrapolated to the range of intermediate energy < 4 GeV) to use it 
as a basic process during the generation of hadron-hadron collisions. After hadronization 
the newly formed secondary hadrons are allowed to rescatter. To determine the interac¬ 
tion between hadrons, the experimental total, elastic, and annihilation cross sections have 
been used. Isotopic invariance and the additive quark model relations were used to avoid 
data dehciency. The resonance cross sections were assumed to be identical to the stable 
particle cross sections with the same quark content. At low energy the QGSM reduces to 
a standard cascade model without mean field effects. The model yields a generally good 
overall £t to most experimental data [24,25]. Particularly the model describes well the 
average kinamatical characteristics and distributions of pions in Mg-Mg interactions [26]. 

We have generated Mg-Mg interactions using Monte-Garlo generator GOLLI, based 
on the QGSM. The events had been traced through the detector and trigger filter. 

In the generator GOLLI there are two possibilities to generate events: 1) at not fixed 
impact parameter b and 2) at fixed b. The events have been generated for b (0< b <6 
fm). The number of simulated events is ~ 4000. From the impact parameter distribution 
we obtained the mean value of < 6 >= 1.34 ± 0.02 fm. For the obtained value of < & >, 
we have generated a total sample of 6200 events. The two regimes are consistent and it 
seems, that in our experiment the value of b=1.34 fm for Mg-Mg is most probable. 

The experimental results have been compared with the predictions of the QGSM for 
b=1.34 fm and satisfactory agreement between the experimental data and the model have 
been found. In Figs 2, 3 and 4 the xp, and - distributions of 7i~ mesons from the 
QGSM calculations are presented together with the experimental ones. One can see, that 
the QGSM reproduces these distributions. The QGSM also reproduces the pp and cosQ 
distributions (Figs.5 and 6). The QGSM data show the similar characteristics in the 
different regions of ( as experimental ones: sharply anisotropic angular distributions in 
the region < (C+ and the almost flat distribution in the region <('■'■>(('+. The slopes of 
Pp - distributions differ greatly in different regions of O'"- The average values < > in 

these two regions also differ: < pp >= (0.029 ± 0.003) (GeV/c)^ in the region Q*" > C"*"; 
< Pp >= (0.109 ± 0.009) (GeV/c)^ in O’" < 0*“- The average values of < > in the 

different regions of ( from the experimental and QGSM data are in a good agreement 
between each other. The distributions obtained by the QGSM in the region > (C+ have 
been fitted by the expressions (8), (9), (10). The results of the fit are given in Table 1 
and Figs. 5 -r 7. In the region C’*' < the pp and the C"'' - distributions have been 
fitted by the formulae (11) and (12), respectively. The results of the fit are given in Table 
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2 and Figs.5 and 7(by the Boltzmann distribution). One can see from the Table 1, that 
the values of the T extracted from the experimental and QGSM data coincide within 
the errors, only for the (1/vr) ■ dN/d(~^ distribution the QGSM slightly overestimates the 
experimental result (for the B-E distribution). 

In our previous article [27] the 7r“ mesons temperature in Mg-Mg collisions was es¬ 
timated by means of inclusive kinetic energy and transverse momentum spectra for tt” 
mesons: in rapidity interval 0.5 -i- 2.1, which corresponds to the pionization region for 7i~ 
mesons and with the c.m.s. angles 90° ± 10°. The pion spectra have been htted by a sum 
of two exponentials, or two temperatures Ti and T 2 . Ti = 55±1 MeV, T2=113±2 MeV. 
This explained by two mechanisms of pion production: directly (T 2 ) and via A resonance 
decay (Ti). The relative yield of T 2 is ~ 22%. One can see, that the light front analysis 
gives temperature which is weighted average value of Ti and T 2 . The temperatures of 
pions have been extracted in the GSI experiments (FOPI, FRS, KAON and TAPS- Gol- 
laborations, see, e.g. [28,29]) similarly. At FOPI Gollaboration [28] it has been obtained 
that the 7r“ spectra from Ni-Ni collisions require the sum of two exponential functions 
with independent yields and slope parameters T; and describing mostly the low and 
the high momentum part of the spectrum, respectively: at E=1.06 A GeV T; = 55 ± 3 
MeV, Th = 93±5 MeV; at E=1.45 A GeV T; = 56±3 MeV, = 100±5 MeV; at E=1.93 
A GeV T; = 61 ± 3 MeV, = 115 ± 6 MeV. At FRS Gollaboration it has been obtained 
that the T for 7i~ mesons in Ne-NaF collisions range from 78±2 MeV to 96±3 MeV 
for projectile energies from 1.34 to 1.94 A GeV. At TAPS Gollaboration for mesons: 
T = 83 ± 3 MeV in G-G interactions at incident energy of E=2 A GeV; T = 70 ± 1 MeV 
in Ar-Ga interactions at incident energy of E=1.5 A GeV; At KAON Gollaboration the 
value of T for 7i~^ mesons ranges from 71±3 MeV (at energy E=1 A GeV) to 95±3 MeV 
(at energy E=1.8 A GeV). 

It should be noted that the extraction procedures of T in the light front variables 
and in the GSI experiments are different. It seems to be interesting in this connection to 
perform the light front analysis of the GSI data. 


5. CONCLUSIONS 


The analysis of tt” - mesons from Mg-Mg collisions in the light front variables has been 
carried out. In some region of phase space of 7r~ mesons the thermal equilibration seems 
to be reached, which is characterized by the temperature T = 75 ± 3 MeV. The variables 



used can serve as a possible convenient tool to study hadron-production processes in 
hadron-hadron, nucleus-nucleus and e~ - interactions. 

A remark on the nature of maxima in -distributions is in order. Recently ALEPH 
Collaboration observed the maxima in the ^ - distributions (.^ = —In p/pmax) [30] of 
secondary hadrons in e~ collisions, which coincide to high precision with the predictions 
of perturbative QCD (see., e.g. [31,32]). The accuracy of coincidence increases when the 
next to leading order corrections are taken into account. So the shapes of - distributions 
are related to the details of the underlying dynamics. Similarly, it seems that the maxima 
in -distributions reflect the dynamics of the processes considered. 


The authors express their deep gratitude to J.-P.Alard, N.Amaglobeli, D.Ferenc, 
Sh.Esakia, S.Khorozov, G.Kuratashvili, J.Lukstins, J.-F.Mathiot, Z.Menteshashvili, G.Paic, 
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FIGURE CAPTIONS 


Fig.l Experimental set-up. The trigger and the trigger distances are not to scale. 

Fig. 2 The xp - distribution of vr mesons, o - experimental data, * - QGSM data. 

Fig.3 The - distribution of tt mesons, o - experimental data, * - QGSM data. 

The curve is a result of polinomial approximation of the experimental data. 

Fig.4 The distribution of tt mesons, o - experimental data, * - QGSM data. The 
curve - result of polinomial approximation of the experimental data. 

Fig.5 The pp distribution of tt mesons, o - experimental data for > (^+ ((^+=2.0), 
o - the QGSM data for > ^+, A - experimental data for < (^+, □ - the QGSM 
data for < (^+. The dashed lines - £t of the experimental data by the Boltzmann 
distribution ; the solid lines - fit of the QGSM data by the Boltzmann distribution. 

Fig.6 The cos© distribution of tt mesons, o - experimental data for > (■+ 
(C+=2.0), o - the QGSM data for A - experimental data for Q*" < □ “ fhe 

QGSM data for < (^+. The dashed lines - £t of the experimental data by the Boltzmann 
distribution; the solid lines - £t of the QGSM data by the Boltzmann distribution. 

Fig.7 The (I/tt) ■ dN/d(~^ distribution of n mesons, o - experimental data, the solid 
line - £t of the data in the region ('+>('+ by the Boltzman distribution, the dashed line 
- £t of the data in the region < C+ by the formula (1 — A - the QGSM data. 

TABLE CAPTIONS 


Table 1. The results of the £t of the distributions of tt 

C+ > C+. 


mesons in the region 


Table 2. The results of the fit of the distributions of tt mesons in the region < (^+ 
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Table 1. The results of the ht of the distributions of tt 

c+ > c+. 


mesons in the region 


Type 

of 

events 

Number 

of 

events 

T (MeV) 

dN/dPf 

dN/ deosQ 

I/tt ■ dN/d(~^ 

Boltz 

B-E 

Boltz 

B-E 

Boltz 

B-E 

exp. 

6239 

76 ± 2 

76 ± 3 

75 ± 3 

78 ± 4 

75 ± 3 

75 ± 4 

QGSM 

6212 

77 ± 2 

80 ± 3 

68 ± 5 

69 ± 5 

82 ± 4 

89 ± 3 
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Table 2. The results of the ht of the distributions of tt mesons in the region < (^+ 



dN/ dp^ 

I/tt * dN/d(~^ 

a 

A 

{GeV/c)-"^ 

^2 

{GeV/c)-^ 

n 

exp. 

0.85 ± 0.03 

12.0 ± 0.4 

4.8 ± 0.3 

4.30 ± 0.06 

QGSM 

0.90 ± 0.05 

11.4 ± 0.6 

5.2 ± 0.9 

4.17 ± 0.11 
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